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Abstract 

The universal law for percolation thresholds proposed by Galam 
and Mauger (GM) is found to apply also to dynamical situations. 
This law depends solely on two variables, the space dimension d and 
a coordinance number q. For regular lattices, q reduces to the usual 
coordination number while for anisotropic lattices it is an effective co- 
ordination number. For dynamical percolation we conjecture that the 
law is still valid if we use the number q2 of second nearest neighbors in- 
stead of q. This conjecture is checked for the dynamic epidemic model 
which considers the percolation phenomenon in a mobile disordered 
system. The agreement is good. 
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1. Introduction 

Recently, a universal power law for both site and bond percolation thresh- 
olds was postulated by Galam and Mauger (GM) (IJ, 0. The GM formula 
is 

p° M = Po[(d - l)(q - lT a d b (1) 

with d being the space dimension and q a coordinance variable. For regular 
lattices, q is the lattice coordination number. The exponent b is either equal 
to a for bond percolation or to for site percolation. Three classes char- 
acterized by three different sets of parameters {a;p } were found. The first 
class includes two-dimensional triangle, square and honeycomb lattices with 
{a = 0.3601; p = 0.8889} for site percolation and {a = 0.6897; p = 0.6558} 
for bond percolation. Two-dimensional Kagome and all (hyper-) cubic lat- 
tices in 3 < d < 6 constitute the second class with {a = 0.6160; po = 1.2868} 
and {a = 0.9346; p = 0.7541} for site and bond respectively. The third class 
corresponds to high dimensions (d > 6). There, b = 2a — 1 and p$ = 2 cl ~ 1 
with a = 0.8800 for site and a = 0.3685 for bond percolation. For hypercubic 
lattices (q = 2d) , the asymptotic limit leading term is identical to the Cayley 
tree threshold. For both site and bond percolation it is, 

Pc=~ (2) 

z 

where z is the branching rate of the tree, i.e. z = q — 1. 

In high dimensions, it has been reported that more GM classes should 
be taken into account ||. Extension to anisotropic and aperiodic lattices 
has also been also reported ||. The formula (1) remains valid. However to 
preserve its high accuracy q should be replaced by an effective non-integer 
value q e ff which is different from a simple arithmetic average. 

In this letter, making a simple conjecture we extend the validity of the 
GM formula (Eq.(l)) to dynamical situations like epidemics or contagion 
models. The results obtained fit well to available data. 

2. Epidemic models 

The epidemic model |J considers the Eden growth |J of a phase in a 



3 



disordered medium with static impurities. In this model, a fraction x of par- 
ticles are randomly dispersed on a lattice. These particles act as hindrances 
for the random growth of a spreading phase. For x < x c , the cluster is grow- 
ing forever while for x > x c the growth of the cluster is stopped after a finite 
number of timesteps. This unblocked-blocked growth transition is closely 
related to a site percolation phenomenon since x c = 1 — p c . 

Recently, the dynamical epidemic model |7], U has been introduced in 
order to provide a phenomenological basis for the aggregation of mesoscopic 
impurities along crystal growth surfaces and interfaces ||. In the dynamic 
epidemic model, the particles are mobile hindrances for the random growth 
of a phase. When the growth front reaches a particle, the latter is supposed 
to move to an empty nearest neighboring site in order to minimize its contact 
with the spreading phase if such a move is possible. This rule is illustrated 
in Figure 1. This repulsive dynamical interaction leads to an aggregation 
phenomenon along the growth interface and leads further to some reorgani- 
zation of the mobile medium |§ . The threshold p c for dynamical situations is 
quite different from the static medium. The formation of growth instabilities 
along the propagating front has been recently reported to be the source 
of the organization process. 

The dynamic epidemic model has been numerically studied at d = 2 
[0, |8| and d = 3 [|10j. Moreover, an exact solution has been obtained on tree- 
like structures which are illustrated in Figure 2. The percolation thresholds 
for either static or dynamic epidemics are listed in Table I as well as the 
thresholds predicted by the GM formula. 

3. The conjecture 

Let us consider first a square lattice. As described above and illustrated 
in Figure 1, when a particle is reached by the interface, the particle tries to 
move to another neighboring empty site as in avalanche processes in sand- 
piles. If the neighborhood of the particle is completely occupied, the particle 
remains static and the cluster growth is pinned there. Intuitively, one can 
consider that the unblocked-blocked growth transition depends at least on 
the occupation of sites in the front neighborhood including the q2 second 
neighboring sites (neighbors of neighbors). Not to be mistaken with next 
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nearest neigbhors. The extended neighborhood of a site is shown in Figure 
3. It is worth noticing that on regular lattices 



q 2 = dq. (3) 

Taking q e ff = q+q 2 = 12, as for percolation with next nearest interactions 
the GM law (Eq. (1), site first class) does not hold and does not reproduce 
the dynamical threshold, for instance p c = 0.44 for squares. Nevertheless, 
one can notice that taking q e ff = q 2 = 8 provides the value 0.4411 for p c 
which is "exact" within the two digits of available data (see Table I). 

Using the same substitution, the universal law holds true also for the 
cubic (d = 3) case for which q 2 = 18 and p c = 015. It yields (for site 
second class) p c = 0.1466 which is again "exact". These results ground the 
following conjecture: Percolation properties of a dynamic epidemic process 
are identical to the associated static ones with substituting q 2 to q. 

We are now in a position to check the above conjecture in the case of 
the Cayley tree for which the thresholds are known exactly. The relevant 
parameter is the branching rate z instead of q. The parameter z can be 
considered as the number of nearest neighboring sites q towards the tree 
extremities. In fact, there is no distinction between directed percolation and 
simple percolation on a tree-like structure. The exact solution of the dynamic 
epidemic model on a Cayley tree |l l| is 



Pc = l/z 2 . (4) 

The number q 2 of second neighboring sites towards the leaves is z 2 . Thus 
from Eq.(2), this theoretical result corroborates the above conjecture that q 
should be replaced by q 2 in the GM universal law. 

On the trees decorated with loops like those illustrated in Figure 2b and 
2c the sites are non-equivalent. An effective number q 2 of second sites to- 
wards the leaves can be calculated by averaging q 2 on the various sites of 
the tree structure. Effective values of the number q 2 for trees decorated with 
loops are given in Table I. Again, considering q 2 instead of q in Eq.(2) pro- 
vides a good value for p c in the case of trees decorated with loops (see Table 
I). However, the agreement is not as good as for the square and cubic lat- 
tices. This discreapancy is to be put in parallel to the finding that for 
anisotropic lattices an effective number of neigbhors should be used instead 
of the arithmetic average. 
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4. Conclusion 

In this letter we have proposed a mapping of the percolation threshold 
for a dynamical situation to a static one by making a simple conjecture with 
respect to the relevant neigbhoring sites. The conjecture states that the 
dynamical percolation case reduces to the static percolation one considering 
that the percolation threshold should include the effect of the second nearest 
neigbhors only. This conjecture was then checked using the GM universal 
law for percolation thresholds. The results are very convincing. It was also 
found to be satisfied in the Cayley tree case. More results on a larger lattice 
spectrum would allow a more definite check of our conjecture. 
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Lattice 


Ref. 


p c static 


p^ M static 


p c dynamic 


pGM dy nam ic 


1 


92 


square (d = 2) 




0.5928 


0.5985 


0.44 


0.4411 


4 


8 


cubic (d = 3) 





0.3116 


0.3115 


0.15 


0.1466 


6 


18 


Tree 


Ref. 


p c static 


l/q 


Pc dynamic 


l/?2 


9 


92 


Cayley tree 


[ 


E3 


l/z 


l/z 


l/z 2 


1/^ 


z 




tree (A loops) 


1 




1/2 


1/2 


0.151 


0.1667 


2 


18/3 


tree (□ loops) 


I 





0.597 


0.625 


0.269 


0.2778 


8/5 


18/5 



Table I — The thresholds p c for both static and dynamic epidemic models 
compared to the associated p^ M thresholds from Eq.(l). All lattices studied 
up to now to our knowledge are listed with pertinent references. The number 
of nearest neighboring sites q as well as the second neighborhood q^ are also 
given. Trees are as those shown in Figure 2. 
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Figure captions 



Figure 1 — Illustration of one growth step of the dynamic epidemic model. 
The growing phase is drawn in white, the mobile particles (hindrances) are 
drawn in black and the empty sites are shown in grey: (a) one empty site 
denoted by a cross in contact with the spreading cluster is selected at random; 
(b) the growth takes place there and the particle touched by the new unit 
jumps towards a neighboring site in order to reduce its contact with the 
cluster. 

Figure 2 — Three different trees with branching rate z = 2 as discussed in 
this paper: (a) Cayley tree, (b) tree decorated with triangular loops, and (c) 
tree decorated with square loops. 

Figure 3 — Nearest neighbors and second neighbors on a square lattice. 
The qi = 8 second neighbors which play the fundamental role in the dynamic 
epidemic model are denoted in grey. 
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